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The mass shift Arrip and width broadening AFp of p mesons produced in heavy ion 
collisions is estimated using a general formula which relates the in-medium mass shift of 
a particle to the real part of the forward scattering amplitude Re/(i?) of this particle on 
constituents of the medium and AF to the corresponding cross section. It is found that 
in high energy {E/A > 100 GeV) heavy ion collisions the p width broadening is large and 
p (or uj) peak could hardly be observed in e~^e~{p^p~) effective mass distributions. In 
low energy collisions {E/A a few GeV) a broad (a few hundred MeV) enhancement is 
expected at the position of the p peak. 

1. INTRODUCTION 

The problem of how the properties of hadrons change in hadronic or nuclear matter 
in comparison to their free values has attracted a lot of attention. Among these proper- 
ties of immediate interest are the in-medium particle's mass shift and width broadening. 
Different models as well as "model independent" approaches were used to calculate these 
effects both at finite temperature and finite density (for a review, see e.g. [||). It is clear 
on physical grounds that the in-medium mass shift and width broadening of a particle 
are only due to its interaction with the constituents of the medium. Thus one can use 
phenomenological information on this interaction to calculate the mass shift and width 
broadening. In a recent paper [Q] two of us argued that the mass shift of a particle in 
medium can be related to the forward scattering amplitude f{E) of this particle on the 
constituents of the medium (in medium rest frame) 

Am{E) = -27r£Re/(E) . (1) 

Here m is the vacuum mass of the particle, E is its energy in the rest frame of the con- 
stituent particle, and p is the density of consituents. The normalization of the amplitude 
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corresponds to the standard form of the optical theorem 



ka = 47rIm/(E) , 



(2) 



where k is the particle momentum. The width broadening is given by 



AT{E) 



P 



ka{E) . 



(3) 



The applicability criteria of Eqs. (1) and (3) were discussed in [2]. In short, they are: 

1) The particle wave length A must be much less than the mean distance between 
medium constituents d, \ = k^^ ^ d. This means that the particle momentum k must 
be larger than a few hundred MeV; 

2) The particle formation length If ~ [E / m) / mchar {f^char ~ "^p) niust be less than 
the nucleus radius R] 

3) Re/(i?) which enters Eq.(l) must satisfy the inequahty | Re/ |< d; 

4) The scattering takes place mostly at small angles, ^ -C 1. Only in this case the 
optical analogy on which Eqs. (1) and (3) are based is correct. 

Eqs. (1) and (3) are correct also in cases when the medium constituents have some 
momentum distributions, e.g. the Fermi distributions for nucleons, or represent a gas at a 
finite temperature. In such cases, in the right-hand sides of Eqs. (1) and (3) an averaging 
over the momentum distribution of constituents must be performed. Eqs. (1) and (3) were 
derived in [2] basing on simple quantum mechanical arguments and the optical analogy. 
This approach allowes one to formulate explicitly the applicability conditions presented 
above. When the medium is a gas in thermal equilibrium the equivalents of Eqs.(l) and 
(3) can be derived basing on thermal field theory. References give few examples and 
the reference gives a relativistic field-theoretic derivation. 

In most of the papers on the in-medium hadron mass shifts the hadrons were considered 
at rest. As seen from Eq. (1) this restriction is not necessary theoretically. It is desirable 
to have theoretical predictions in a broad energy interval, since it extends possibilities of 
experimental investigations. As discussed in for the cases of p or vr-mesons embedded 
in nuclear matter the energy dependence of the mass shifts is rather significant at low 
energies, i.e. in the resonance region. 

We estimate the p meson mass shift and width broadening in the case of p-mesons 
produced in heavy ion collisions. The most interesting case is the case of which can be 
observed through the decay p° e~^e~ {p~^p~). We will assume that p-mesons are formed 
at the latest stage of the evolution of hadronic matter created in the course of a heavy ion 
collision when the matter can be considered as a noninteracting gas of pions and nucleons. 
(We will neglect the admixture of kaons and hyperons, which is known to be small 0, as 
well as heavy resonances.) This stage is formed when, during the evolution of the matter 
after collision, the total density of nucleons and pions becomes of order of the normal 
nucleon density in nuclei. The description of nuclear matter as a noninteracting gas of 
nucleons and pions of course cannot be considered as a very good one. So, it is clear from 
the beginning that our results may be only semiquantitative. The main ingredients of our 
calculation are the pvr and pN forward scattering amplitudes and total cross sections as 
well as the values of nucleon and pion densities. 
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We consider here central heavy ion colhsions and assume that the nucleon and pion 
momentum distributions in the gas are just the momentum distributions measured exper- 
imentally in such collisions. (Another model, in which nucleons and pions are assumed to 
be in the state of ideal gas at fixed temperature and chemical potential, will be considered 
in a subsequent publication [8].) 



2. CALCULATION OF pvr AND pN CROSS SECTIONS AND FORWARD 
SCATTERING AMPLITUDES. 

Let us first focus on the amplitudes and cross sections. To determine these quantities we 
use the following procedure. In the low energy region we saturate the cross sections and 
forward scattering amplitudes by resonance contributions. At high energies we determine 
apN and a p.,, from a^N and a^T, using vector dominance model (VDM). a^N is well known 
experimentally 0, Re/^A? is determined from dispersion relation; ct^tt and Re/^jr can be 
obtained by the Regge approach. Since VDM allows one to find only the cross sections 
of transversally polarized p-mesons, we restrict ourselves to this case. As was shown in 
[2], dX Ep >j 2 GeV Am and AF for longitudinal p-mesons in nuclear matter are much 
smaller than for transverse ones. At zero p-meson energy Am and AF for transverse and 
longitudinal p-mesons are evidently equal. In the case of scattering on low temperature 
gas they are comparable [10]. Therefore, for unpolarized p-meson our results should be 
multiplied by a factor ranging from 1 to 3/2. 

To estimate Re/p7r(s) at low energy we write in the center of mass (cm.) frame 

Re/,.(.) = - Y. F.F.^ ,^«'"'^r"?',, . (4) 

R ^Qcm Ws - niRY + F|j/4 

where -y/i is the total cm. energy, m/j and T r are the mass and the total width of the 
resonance, Br is the branching ratio of its decay into Trp and Qcm is the center of mass 
momentum, 

qcm = {mp + m^y][s - {mp-m^y]. (5) 

Z^y S 

The factor Fg is the spin factor, Fi is the isospin factor. The latter is equal to 1/3 = 
(1/2) X (2/3) for Ir = 1. The first factor here refiects the fact that we are interested 
in the p^vr scattering and only one of the two decay channels of a Ir = 1 resonance 
can contribute here: — ^ P°7i"^, but not p^7r°. The second factor corresponds to the 
assumption, that all three pion isospin states are equally populated in the gas. Similarly, 
for Ir = the isospin factor is 1/9 = (1/3) x (1/3). We take into account the following 
resonances [9]: ai(1260), 7r(1300), 02(1320) and ti;(1420). The nearest resonance under 
the threshold, to'(782), does not give a considerable contribution due to its narrow width. 
For the spin factor we take Fg = 1, 1, 2, 1 correspondingly for the above mentioned 
resonances. (These factors take into account that we consider only transvers p). The 
amplitude in the pion rest frame is obtained from Eq. (4) multiplying it by the rescaling 
factor kp/qcm = Vs/mT^, where kp = J E"^ — rrip is the p momentum in the pion rest frame. 
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For Gfyjr we use the standard resonance formula 



' fm 



TT 



R 



(6) 



As is known, the pion scattering amphtude on any hadronic target vanishes at zero 
pion energy in the target rest frame in the hmit of massless pions (Adler theorem). In the 
framework of effective Lagrangians this can be achieved if the pion field enters through 
its derivative, dip/dx,j,. We assume that in the interaction describing the fm scattering 
through the ai resonance, dip/dx^ is multiplied by the p-meson field strength tensor F^^ 
and the field ai^. This results in appearance of an additional factor in Oi contributions 
to Re/pTT and cXpTr in Eqs. (4) and (6) 



(7) 



which is normalized to 1 at s = m^^. At s > m^^ this factor was not taken into account. 
Similar factors were also introduced for contributions of other resonances. 

In the high energy region we assume that the Regge approach is valid for jti scattering 
and apply VDM to relate pTi and 77r amplitudes. As is well known, the Regge pole 
contributions to the forward scattering amplitude normalized according to Eq. (2) have 
the form 



Atts 



k ^ 1 + e 



sm irai 



(8) 



where ctj is the intercept of i-th Regge pole trajectory, is its residue, k is the projectile 
momentum in the target rest frame. As follows from Eqs. (2) and (8), 



a{s) = J2 riS 



(9) 



Ti \ ^ V- l+COSTTQij 

Re/(s) = - — E — -■ '^'■^ 



sm ttq;,- 



(10) 



For (T^^ only the P (pomeron) and P' poles contribute [11,12]. The residues of P and 
P' poles in 77r scattering were found by Boreskov, Kaidalov and Ponomarev (BKP) [12] 
using the Regge pole factorization theorem and the data on 7p, irp and pp scattering. 
Taking the BKP values for the P and P' residues we have 



(77r-y(s) = 7A8a 



^SoJ \SoJ 



(11) 



where ap = 1.0808, ap> = 0.5475, a = 1/137, sq = 1 GcV^ and a in Eq. (11) is given in 
millibarns, if s is in GcV^. For the P and P' intercepts we take the Donnachic-Landshoff 
values [13]. Since in their fit to the data BKP assumed ap = 1 and ap' = 1/2, the 
values of residues in Eq. (11) are shghtly changed in comparison with [12] in order to get 
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the same Gt^^ at s = 9 GeV^. From Eqs. (10) and (11) the real part of the forward 77r 
scattering amphtude can be found: 




Re^,(s)^ rest frame = -—7.48a <^ -0.106 - + 0.752 - \ , (12 



where the momentum k is in GeV and Re/ is in mb-GeV. 

In VDM (Tp7r(s) is related to <7^tt{s) by (see [2]) 




cTp^{s) = -^{l + ^\ a^^is), (13) 



where gp/4:7T = 2.54, Qp/g"^ = 1/8 and the <^-meson contribution is neglected. A similar 
relation holds also for Re/p,r- Unlike in [2], we prefer here to use direct Regge formulae 
for Re/ at high energies instead of finding it from a through the dispersion relation, since 
in the latter approach the results are sensitive to the low energy domain which is more 
uncertain. 

The results of the calculations of ap-„ and Re/p7r as functions of p-meson energy in the 
pion rest frame are presented in Fig.l. As seen from Fig.l, the matching of low and high 
energy curves is satisfactory. 

For the amplitude RefpN at laboratory frame energies of p above 2 GeV we use the 
results of Ref. obtained using dispersion relation, VDM and experimental data on a^^. 
At lower Ep we again use the resonance approximation 

(the factor (1/4) appears because we consider only transverse p). The isospin factors Fi 
are 1/3 and 2/3 correspondingly for and A resonances. We take 10 well established 
A^ and A resonances [9] with significant branchings into pN and with masses above the 
pN threshold and below 2200 MeV. The set of baryonic resonances taken into account is 
close to the set used in [14] . The main difference in comparison with [14] is that effective 
widths = [qcm/ <lcm,Rf'''^^ were introduced only for the resonances close to the pN 
threshold {qcm,R is the value of qcm at the resonance). At qcm > QcmR we put Fg^ = F^^. 
Apart from these resonances, in the calculation of Re/pA? two resonances with masses 
below the pN threshold were also accounted for: A(1238) and A^(1500). It was assumed 
that VDM relates TpN and T^n of these resonances in the following way. Since both 
resonances are close to the pN threshold, we can write for each of them: Fp^r = qcmlpN, 
T-yAT = kcml'yN, whcrc qcm and kcm are the pN and 7A^ cm. momenta, respectively. Then 
we assume that 7pAr and j^n are related by the VDM formula 



7^Ar = 47ra— 1 + -| 7pAr . (15) 




The values of 7^Ar can be found from the values of ct^at at the resonance peaks. The 
contributions of A(1238) and A^(1500) to RefpN are essential at low energies of p mesons: 
they comprise about (—1) — (—0.5) fm at ii^p = 1 — 2 GeV in the nucleon rest frame. 
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The results for and HefpN in the rest frame of the nucleon together with the curve 
obtained in Ref.0] for high energies and the matching curve are shown in Fig.2. As can 
be seen, the matching of low energy and high energy curves is good. 



3. DETERMINATION OF THE MASS SHIFT AND WIDTH BROADEN- 
ING OF THE p-MESON PRODUCED IN HEAVY ION COLLISION FROM 
THE DATA ON NUCLEON AND PION DISTRIBUTIONS 

As was mentioned above, in heavy ion collisions we will consider only nucleons and 
pions as constituents of the medium. Therefore, in this case Eqs. (1) and (3) take the 
form 

Amp{E) = -— {p^Refpr^iE) + p,Re/,,(E)} , (16) 
m 

ATp{E) = ^ {pN apN{E) + p^ap^{E)} , (17) 

where p^r and p^^ are the nucleon and pion densities in the medium formed at the latest 
stage of evolution of hadronic matter produced in heavy ion collisions. 

We will restrict ourselves to consideration of head-on (central) collision with small values 
of impact parameter when the number of participants - the nucleons which experience 
considerable momentum transfer - is close to the total number of colliding nucleons. 

Experimental data shows that the nucleons and pions produced in heavy ion collisions 
cannot be considered as thermal gas even at the final stage of evolution of hadronic matter 
created in the collisions. In order to demonstate this let us discuss separately the cases of 
high energy {E/Ar^ 100 GeV) and low energy {E/Ar^a, few GeV) heavy ion collisions. 
In case of high energy collisions the longitudinal (relative to the beam) and transverse 
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momenta of nucleons and pions are very different. In the experiment on S + S collisions 
at E/A = 200 GeV [15] it was found (in the centre of mass frame) that (pl^) — 3-3 GeV, 
(Ptn) = 0.61 GeV and (p^™) ~ 0.70 GeV, (pt-k) ~ 0.36 GeV. In other experiments on high 
energy heavy ion coUisions (see e.g. [16,17]) the situation is qualitatively similar. This 
means that in this case one can by no means speak about thermal gas of final particles 
and their momentum distributions must be taken from experiment. 

The data on low energy heavy ion collisions also indicate that pions and nucleons cannot 
be described as gases in thermal equilibrium. The angular distributions of pions produced 
in Ni + Ni collisions a,t E/A — (1 — 2) GeV show a considerable anisotropy[18]. If the 
pion angular distribution in the centre of mass frame is approximated by 1 + a cos^ 9, then 
it follows from the data that a ^ 1.3. Unfortunately, there is not enough information 
in the data about the nuclcon angular and momentum distributions. We have checked 
the hypothesis of thermal equilibrium by assuming that the probability of production 
of a given number of particles is proportional to the statistical weight of the final state 
(Fermi-Pomeranchuk approach [19,20]). It is evident that this hypothesis is even more 
general than the hypothesis of thermal equilibrium. In this approach the pion/nucleon 
ratio -Rtt = ^t^/N in central collisions can be predicted in terms of the main ingredient 
of the method - the volume per particle at the final stage of evolution (and, of course, 
the initial energy). The calculations show that the data [18] on the energy dependence of 
the ratio Rt^ are well described by the statistical model, but in order to get the absolute 
values of R-,^ in Ni + Ni as well as in Au + Au [21] collisions, it is necessary to put the 
volume per nuclcon very small, about 5 times smaller, than the one in a normal nucleus, 
which is unacceptable. 

Therefore, the only way to perform the averaging over momentum distributions of pions 
and nucleons seems to take the latter from experimental data on heavy ion collisions. 

In calculation of the p-meson mass shift and width broadening the averaging must 
be performed over the p-mcson direction of flight relative to nucleons and pions. Such 
calculation can be done only for real experimental conditions. For this reason we restrict 
ourselves to crude estimates. 

Consider first the case of high energies and take as an example the experiment [15] 
on central collisions, where the ratio of pion to nucleon multiplicities was found to be 
-Rjr = 5.3 . Suppose that in this experiment the p-meson is produced with longitudinal 
and transverse momenta in the lab. frame k^p = 3 GeV, = 0.5 GeV. We choose these 
values as typical for such experiment. At such values of p momenta the formation time of 
p-meson is close to the mean formation time of pions produced in heavy ion collisions. So, 
a necessary condition of our approach is fulfilled. Since the mean momenta of nucleons 
and pions in the experiment [15] are known (they were presented above), it is possible 
using the curves of Figs. 1 and 2 to calculate the mean values of Refpj^, Refp^^, apN and 
cTpTT in pA^ and pvr scattering. The results are (in the lab. frame): 





(18) 



(cTpjv) ~ 43 mb , 



(cTpTr) ~ 20 m6 . 



(19) 



The small value of (Rc/p^) comes from the compensation of positive and negative Re/p^r 
at low and high energies (Fig. 1), i.e. from the scattering of p on pions moving in the same 
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direction as p (comovers), or in the opposite one. Because of this compensation (RefpT^) 
is badly determined, but since it is small this does not affect the final result. Using Eqs. 
(16) and (17) we can now find the mass shift and width broadening of p- meson. For the 
nucleon and pion densities we take 

_N_ N _ 1 

n n 1 , . 

P^'^V^ Nv^ + nv^ ^ ^%^(l + i4^) ' ^ ^ 

where N and n are the numbers of nucleons and pions at the final stage of evolution, 
Rtt — n/N, V is the volume of system at this stage. It is assumed that at this stage any 
participant (nucleon or pion) occupies a definite volume or v^^. We can write 

'^ = '-TTr;p^ p^ = p-''-iVr7p^ ^''^ 

where p^vr — and [3 = v-^/v^. For numerical estimates we take p% = 0.3 fm~^, 

almost two times the standard nucleous density. This number is probably one of the most 
uncertain ingredients of our calculations. Substitution of Eqs. (18), (19) and (22) into 
Eqs. (16) and (17) gives (for the experimental value of Rt^ — 5.3 and P — 1) 

Arup = 18 - 2 = 16 MeV , (23) 

ATp ^ 150 + 400 = 550 MeV . (24) 

The first numbers in Eqs. (23) and (24) refer to the contributions from pN and second 
ones from pn scattering. Because p-meson width broadening appears to be very large, 
the basic condition of our approach, ATp <^ rrip is fullfilled badly. The other applicability 
condition of the method, | Re/ |< d, is also not well satisfied, since in this case d = 0.9 
fm. For these reasons the values of ArUp in Eq. (23) and of ATp in Eq. (24) may be 
considered only as estimates. 

The main conclusion from Eqs. (23) and (24) is that in the case of p-mesons produced 
in high energy heavy ion collisions with the longitudinal and transverse momenta chosen 
above, the mass shift is small, but the width broadening is large and hardly a p (or lu) 
peak would be observed in e+e" or p"^p~ effective mass distributions. Let us estimate 
how sensitive are the results to a variation of /c^ and kr- It can be easily seen that the 
mass shift will be small in all cases (say, Arrip < 50 MeV). If we put kxp = instead of 
krp = 0.5 GcV, this will only weakly affect the mean value of CTpi^ and decrease (Tp^^ by 
20%. The latter results in decrease of ATp in Eq. (24) by 80 MeV, i.e. within the limits 
of accuracy. The variation of kpL between 1 GeV and 10 GeV also results in variation of 
the same order (10-20%) in AFp in Eq. (24). 

As was mentioned above, the main uncertainty in our approach comes from the value 
of the nucleon density at the final stage of evolution which we put to be p% = 0.3 fm^'^. If 
this density at the time of the p-nieson formation would be, say, two times smaller, then 
we would have AFp ~ 250 MeV and the p-meson could be observed as a broad peak in 
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e"'"e~ or mass spectrum. It should be mentioned, however, that the chosen above 

value of f3 = Vt^/vj^ = 1 is not the only possible, or even the most plausible. If we assume 
that (3 = (r^/rAr)^, where r^r and tn are pion and nucleon electromagnetic radii, r^r = 0.66 
fm, rjv = 0.81 fm, then /3 Ri 0.55. This choice of P increases ATp by a factor of 1.6. 

Our qualitative conclusion is the following. In the central heavy ion collisions at high 
energies, E/A 100 GeV, where a large number of pions per participating nucleon is 
produced, the p- (or u) peak will be (if at all) observed in e^e" or fi^fi^ effective mass 
distributions only as a very broad enhancement. Inspite of the fact that in our approach 
we used a few assumptions, which may raise some doubts (hypothesis of noninteracting 
nucleon and pion matter at the final stage of evolution, the numerical value of nucleon 
density, etc.), we believe that this qualitative conclusion is still intact. 

This conclusion is in qualitative agreement with the measurement of e+e^ pair produc- 
tion in heavy ion collisions [22], where no p peak was found and only a smooth e+e~ mass 
spectrum between and 1 GeV was observed. If, however, such a peak will be observed in 
future experiments, this would indicate that the hadronic (nucleon and pion) density at 
the final stage of evolution, where p-meson is formed, is low, lower even than the normal 
nuclear density. 

Recently preliminary data have been presented [22], where in studying the e+e~ mass 
spectrum in Pb - Au coUisions a,t E/A— 160 GeV it was found that the p-peak is absent 
for A;T(e+e~) < 400 MeV, but reappears as a broad enhancement for A;r(e+e~) > 400 
MeV. We do not see a possibility for this in the framework of our approach in the case of 
central heavy ion collisions. Moreover, we believe that for central collisions the absence 
of p-peak at low /cr and its reappearance at higher will be hard to explain in any 
reasonable model. The only explanation we see for this effect (if it will be confirmed) , is 
that in this experiment the peripheral p-meson production plays an essential role. In this 
case the p- mesons with higher kr have a larger probability to escape the collision region 
and decay as free ones. 

Let us turn now to the case of low energy heavy ion collisions, E/A a. few GeV. 
Consider as an example the case of heavy ion collisions at Ej^^n/A — 3 GeV and production 
of p-meson with the energy = 1.2 GeV in the forward direction. (This energy of p- 
meson was chosen, because our approach works better at higher Ep, and p of this energy 
can be kinematically produced at such heavy ion energy). The number of pions, produced 
in Ekin/A = 3 GeV collisions can be found by extrapolation of the data [18] on Ni + Ni 
collisions. The data shows that i?^ is with a good accuracy linear in y/s/2 — m. We 
find i?7r = 0.48. As follows from the analysis of the data [18] at E^in/A — 1.93 GeV, 
the energies of produced pions are rather small, E^, ~ 200 — 300 MeV . At such low 
energies it is reasonable to suppose for pions < >=< p^^^ >k, {).2 GeV . Assuming 
(this assumption does not influence essentially the final results) that the mean transverse 
momentum of nucleon participants produced in the collision is the same as at high energy, 
< Ptn >— 0.61 GeV [15] (see above), we can construct the momentum distributions of 
nucleons. Then we are in a position to calculate the mean values of RefpN, Refp^^, ap^ 
and cr„7r for this case. The results are: 



(Re/p;v) = -0.54 fm, 

Wpn) =4:6mb, 




(Refp^) = 0.30 fm , 



(25) 

(26) 
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For the p mass shift and width broadening we have (at the same value of p% as before 
and 13 = 1): 

Arup = 37 - 10 = 27 MeV , (27) 

ATp = 250 + 35 = 285 MeV . (28) 

(The first numbers in Eqs. (28) and (29) refer to pN scattering, the second ones to pn). 
The conclusion is that at low energy heavy ion collisions p-peak may be observed in e~^e~ 
or fj,~^fi~ effective mass distribution as a broad enhancement approximately at the position 
of p-mass. 

We are indebted to K. Boreskov, A. Kaidalov, G. Brown and A. Sibirtsev for illumi- 
nating discussions. B.I. thanks A. Smirnitsky and V. Smolyankin for the help in getting 
information about experimental data. 
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